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■ Abstract. The Spitzer Space Telescope and in particular its InfraRed Array 
Camera (IRAC) is an ideal facility to study the distribution of AGB stars in our 

^ , own and other galaxies because of its efficiency in surveying vast areas of the sky 

• and its ability to detect sources with infrared excess. The IRAC colors of AGB 
[ stars, however, are not well known because cool stars have numerous molecular 

absorption features in the spectral region covered by the IRAC photometric 
system. The presence and strength of these features depends on the chemistry 
of the stellar atmosphere and the mass loss rate and can change with time due 
to the star's variability. To characterize the IRAC colors of AGB stars, we 
^ ' are carrying out a Spitzer Guaranteed Time Observation program to observe 

i a sample of AGB stars with IRAC. The results will be made available to the 

■"sj" . community in the form of template magnitudes and colors for each target with 

' the goal of aiding the identification of AGB stars in already available and future 

IRAC surveys. We present here the first results of this project. 

o 

J2 ■ 

P^. 1. Mapping AGB Stars with Spitzer /IRAC 

2 . Asymptotic Giant Branch (AGB) stars, with luminosities exceeding 10^ Lq, 

' are among the brightest stars in the galaxy. They are also among the reddest 

• because of intense mass loss processes (up to 10~^ Mq/jv) responsible for en- 
shrouding AGB stars in dusty envelopes, which are the source of strong infrared 
excess. These characteristics of AGB stars make them very important tools to 

^ I study the structure of our and other galaxies by mappin g their distribu tion at 

■ infrared wavelengths. The Spitzer Space Telescope (Werner et al.|[2003 ) is es- 



pecially suited for this task. The telescope has ver y low background emission, 
which makes the InfraRed Array Camera (IRAC, iFazio et al.l |2004| ) onboard 
Spitzer very sensitive and allows mapping large areas of the sk y in a very short 



time . Several large area IRAC surveys, including GLIMPSE (iBeniamin et al 



2003) covering 220 square degrees of the galactic plane and SAGE (" Meixner et al 



2006 ). mapping a 50 square degrees area centered on the Large Magellanic Cloud, 



are already available. The four IRAC channels operating at 3.6, 4.5, 5.8, and 
8.0 /im allow detection of both photospheric and dust shell emission. The IRAC 
colors of AGB stars, however, are not very well known. 

AGB stars have numerous molecular abso rption features in the spectral 
region covered by the IRAC photometric system ( Waters et al.llT999l ). a s demon- 



strated by observations with the Short Wavelength Spectrometer (SWS. IValentiin et al 



[1996) aboard the Infrared Space Observatory (ISO). Features seen include H2O, 
SiO, CO2, CO, and silicates in stars with atmospheric C/0 ratio <1, while 
C2H2, HCN, CS, C3, and carbonaceous dust (SiC and amorphous carbon) have 
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been found in carbon stars, which have C/0 > 1. The presence and strength 
of these s pectral features depends on the chemistry of the stellar atmosphere 
(jSloan et al . 1_998; Slo aii &: Price 199p.) and the mass loss rate. The features can 
also change with time ( Qnaka et al. 20021 ) because AGB stars are long period 
variables of Mira, semi-regular, or irregular type. 

As a result of these dust and molecular features, the IRAC colors of AGB 
stars can be quite different from the "reddened photospheres" that one would 
expect for mass losing giants of late spectral type. There is a need, in order to 
efficiently identify AGB stars among other red objects, for accurate measure- 
ments of IRAC colors of AGB stars. For this reason, we have started a program 
to observe a sample of nearby AGB stars with IRAC as part of Cycle-3 Guar- 
anteed Time Observations (GTO). The program is in progress, and we present 
here the first results. 



2. IRAC AGB Colors from ISO SWS and Models 



As a first step in deriving IRAC colors for AGB stars, we computed synthetic 
colors for 87 0-ri ch AGB stars, 27 carbon stars, 8 S-sta rs, and 23 M-type super- 
giants (lists from lSloan k Pricelll998l :[sl oan et allll998l l by convolving their ISO 
SWS spectra with the IRAC bandpasses. We obtained the ISO spectra (2 12 in 
total, with some stars observed in multiple epochs) from lSloan et al. ( 20031 ) and 
the IRAC transmission profile from the Spitzer Science Center web sit^B- 

The distribution of the sources on two IRAC color-color diagrams is shown 
in Figure 1: 



1. Sources with larger mass loss, as expected, have a larger infrared excess, 
especially in the [3.6] — [8.0] color that is the most sensitive to the overall 
slope of the spectrum. 

2. IRAC colors are not sensitive to the detailed dust composition because even 
the longest passband, centered around 8.0 fim, only marginally includes 
the 10 fim silicate feature. 

3. Carbon stars tend to have redder [3.6] — [4.5] colors than 0-rich stars 
having the same infrared excess. This is due to very strong C2H2 and 
HCN absorption within the 3.6 fim band in carbon stars and broad CO2 
and SiO absorption features in the 4.5 ^m band in 0-rich stars. 

4. There is a population of carbon stars with very blue (—0.5 magnitudes 
below the model tracks) [4.5] — [5.8] color due to strong C3 absorption at 
5.8 fim. This feature is transient, as some sources observed by ISO at 
multiple epochs have this anomalous color only at some pulsation phases. 

5. S-stars and M-type supergiants tend to have colors very similar to the 
colors of M-type AGB stars. 



^http://ssc. spitzer. caltech.edu/irac/spectraLresponse. html 
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Figure 1. IRAC synthetic color-color diagrams of O-rich AGB stars (cir- 
cles), carbon stars (squares), S-stars (triangles), and M-type supergiants (as- 
terisks) derived from ISO SWS spectra. The tracks are models of circumstellar 
envelopes with silicate (solid lines) and amorphous carbon (dashed lines) dust. 
Thin lines are DUSTY models usi ng a black body as ce ntral source, and the 
thick lines are models computed bv lGroenewege using realistic stellar 

atmosphere models. The two O-rich outliers are two-epoch measurements of 
V354 Lac, an AGB star with a very thick envelope. 
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Figure 1 also shows the tracks of simple radiative transfer models in the 
IRAC colors. The DUSTMI models, which use a cool black body for the spec- 
trum of the central source, cannot fit the colors of sources with little or no 
infrared excess because the models are missing the i nolecular features res ponsi- 
ble for blue IRAC colors. The models computed by iGroenewegen usmg 
realistic spectra for the central AGB star fit the data much better. On the whole, 
the model tracks intercept the synthetic colors for most of the sources with the 
exception of the carbon stars that show strong C3 absorption features. 



3. The Spitzer/IRAC GTO Program 

The nearby AGB stars selected for IRAC observations are ones for which ac- 
curate parallax or interferometric distances and reliable determinations of their 
mass loss rates are available. The observations target 22 0-rich stars, 7 intrinsic 
S-stars, 19 carbon stars, and 4 M-type supergiants. The supergiants have been 
added as a comparison sample of mass- losing stars outside the AGB. The mass 
loss rates of our target stars range from 10~® to 10~^ Mq/jt in each category 
( Guandalini et al. 20061 ). The sample contains Mira, semi-regular, and irregular 



variables. 

The goal of the program is to measure accurate IRAC photometry for all 
stars in the sample in order to validate and cross-calibrate the ISO SWS spectra 
for this class of sources. This is particularly important because the absolute 
photometric calibr ation of IRAC is based mainly on A-type primary calibrators 
( Reach et al.ll2005l ) and because of the uncertainties in the ISO SWS absolute 



calibration resulting from the splicing of the different spectral segments. Each 
star will be observed in two epochs two months apart in order to evaluate the 
change in color at different variability phases. (The timing of the observations 
is constrained by the orbit of the spacecraft). 

Nearby AGB stars are of course among the brightest objects in the infrared 
sky. This makes the IRAC photometry very difficult because all sources will 
be heavily saturated even with the shortest available frame times. To solve 
this problem, we have developed a technique to derive the Vega magnitudes of 
saturated stars by fitting the low level features of the saturated PSF (diffraction 
spikes and P SF wings) with a mo del of the IRAC PSF derived from a sample of 
bright stars ([Marengo et al. 20061 ). Figure 2 shows the 4.5 /um image of SZ Car 



(SRb carbon star), one of the first objects observed in our program, before 
and after PSF subtraction. The stability of the Spitzer optical system allowed 
us to create a very accurate model of the IRAC PSF. The PSF is directly 
normalized to the actual image of Vega, which is one of the stars used in its 
construction. By fitting the observed saturated sources with this PSF, we can 
derive their magnitudes with an accuracy within 1 - 3% independently of the 



standard IRAC flux cahbration (jSchuster. Marengo fc Patten 2006). 



The program is in progress. As of 2006 November, 33 stars have been 
observed in their first epoch. Preliminary results show that the photometry is 
consistent with the synthetic colors derived from ISO SWS. The program will 



^http: //www. pa.uky.edu/~moshe/dusty 
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Figure 2. IRAC 4.5 /um image of SZ Car (SRb carbon star with m[4 5] = 
1.18). Left panel is the mosaic resulting from the co-addition of 5 individual 
2-second frame time exposures. The saturated star occupies almost all the 
IRAC field of view (5x5 arcmin). The right panel shows the same image after 
PSF subtraction. The PSF fitting determines the magnitude of the AGB star 
with an accuracy of 1 - 3% in flux. 



be completed within one year, after which we will release the complete catalog 
to the community. 
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